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Cp  Specific  heat  of  steel 

D  Diameter  of  ingot 

H  Heat  of  fusion  of  steel 

h^  Surface  coefficient  of  heat  transfer,  ingot  to  mold 

hyf  Surface  coefficient  of  heat  transfer ;»  water  side 

X  Mold  wall  thickness 

L  Hold  len^ 

q  Heat  transferred,  Btu/^r«,ft^ 

Q  Heat  transferred,  dtu/h^'* 

ri  radius  of  ingot 

T2  radius  of  mold,  water  side 

t^  average  surface  tei^erature  of  ingot  within  mold 

t2  average  water  temperature  in  mold 

T  mean  temperature  of  ingot  cross-sectiai  at  mold  eacit 

Tp  pouring  temperature 

U  Composite  heat  transfer  coefficient 

V^  velocity  of  ingot 

^  thermal  diffusivity 

b  density  of  steel 

hwt>/j^  Nusselt  number,  water 

^*^//Af  iteynold*s  number,  water 

^Jy\s.  Prandtl  number,  water 
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SUMMABT 

AaalyBis  of  continaous  casting  involyes  many  piiases,  two  of  which 
ara  described  in  tiiis  thesis*  The  first  investigation  consists  of 
detemining  the  temperature  distribution  in  a  moving  ingot  of  metal 
which  is  solidified  and  cooled  sufficiently  in  a  water  cooled  mold 
to  peroit  the  continaous  processing  of  the  resulting  billet*  In 
the  second  phase  there  has  been  establi^ed  a  relationship  for  Buudimim 
water  flow  rates  for  heat  transfer  and  their  association  with  cor- 
responding  coefficients  of  heat  transfer. 

The  analysis  has  ^own  that  the  rate  of  heat  transfer  is  primarily 
a  function  of  the  siirface  coefficient  of  heat  transfer  between  the 
ingot  and  the  mold*  To  increase  this  rate  some  means  must  be  provided 
to  increase  this  coefficient.  Addition  of  oil  along  the  boundary  has 
bacda  suggested  xriiich  would  not  only  increase  this  surface  coefficient 
but  will  also  decrease  the  frictional  drag* 

The  results,  although  only  approxixaate,  are  provided  as  a  guide 
for  future  work* 
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Generally,  the  analytic  mode  of  approach  to  the  solution  of  a 
problem  also  requires  exj^rimentation.  However,  lack  of  facilitieo  and 
time  limitations  precluded  experiawatation  for  verification  of  any 
results  in  this  thesis. 

The  priaary  endeavor  in  this  investigation  is  to  assist  the 
experimenter  in  his  woric  by  forming  limits  of  practicability  and 
providing  raathnatical  solutions  to  in^rtant  x^iases  of  the  problasu 
Results  are  only  approximations,  consistent  with  assumptions  necessary 
to  perform  the  operations,  but  they  can  be  used  as  a  guide  for  future 
work. 

Ti^atment  of  the  entire  problem  of  continuous  casting  is  very 
long  and  complex,  and  to  the  best  1a»3wledge  of  the  writer  no  work  has 
been  done  analytically  along  these  lines,  either  in  full  or  upon  any 
lidiase.  No  information  is  available  for  many  of  the  problems;  however, 
with  work  being  done  at  the  present  tiiae  on  properties  of  metals  at 
temperatures  at  and  near  the  melting  point,  some  hope  is  held  that  in 
the  future  information  will  be  available  concerning  these  properties. 

The  process  involves  pouring  molten  steel  (with  a  minimum  amount 
of  superheat  to  allow  good  fluidity)  into  a  vertical,  water  cooled, 
metal  mold.  At  the  start  of  the  process  the  mold  is  plugged  at  the 
lower  end  until  sufficient  freezing  takes  place  to  allow  adequate 
strength  to  permit  pulling  the  plug  and  continuously  processing  the 
resulting  billet. 

An  analysis  involves  investigation  of  various  phases  which  are 
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omtXined  below  and  discussed  later  In  this  papers 

(a)  Teasperature  distrilmtion  at  any  point  in  the  ingot* 

(b)  Investigation  of  limiting  water  flow  rates,  and 

(o)  Prediction  of  accompanying  coefficients  of  heat  transfer 
between  the  ingot  and  the  mold* 

(d)  Possibility  of  frictional  drag  exerted  on  the  mold  by 
the  ingot* 

(e)  !4axijmim  velocity  of  ingot  travel* 

(f)  Load  capacity  of  the  ingot* 

(g)  Deteiroination  of  ingot  breaking  point* 

(h)  Investigation  of  possible  surface  boiling  on  water  side 

of  the  mold* 
(i)  Investigation  of  surface  tension  and  surface  film  effects* 
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The  greatest  amount  of  vork  on  continuous  casting  of  steel  has 
been  done  by  Edsrard  R,  Williams  and  the  Williams  Engineering  Company, 
His  solutions  to  the  problems  have  been  of  an  experij^ntal  natiire, 
and  to  the  best  knowledge  of  the  writer,  no  attempt  has  been  mad© 
toward  an  analytical  investigation. 

The  WilHaas  process  accomplishes  the  necessary  heat  removal  by 
providing  a  very  thin  mold  wall  of  high  conductivity  material  (cobh 
mercial  drawn  copper  or  brass  tubing)  throng  which  heat  can  rapidly 
pass  from  the  ingot,  and  by  removing  this  heat  extremely  rapidly  from 
the  outside  of  the  mold  wall  by  means  of  water  or  other  cooling  fluid 
passed  at  hi^  velocity  along  the  surface.  By  removing  the  heat 
rapidly  from  the  outer  surface  of  the  iiwld  wall,  a  hi^  temperature 
gradient  is  produced  in  the  mold  wall,  which  permits  fast  removal  of 
heat  from  the  ingot.  It  is  claimed  that  this  instantly  freezes  the 
molten  metal  into  a  skin  which  shrinks  away  from  frictional  resistance 
and  is  sufficiently  strong  to  prevent  fracturing  as  it  is  moved  along 
the  mold  wall. 

The  cooling  water  in  the  mold  passes  throu^  passages  which  have 
a  very  small  clearance,  the  velocity  being  quite  hi^.  It  has  been 
found  that  a  relatively  small  amount  of  water  is  sufficient  to  satis- 
factorily remove  the  heat  from  a  steel  ingot  to  enable  it  to  be  poured 
and  %ri.thdrawn  at  a  speed  of  seven  or  more  feet  per  minute. 

As  the  ingot  emerges  from  the  bottom  it  passes  througti  a  series 
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of  sizing  rolls  that  prevent  bulging  of  the  walls  by  ferrostatic  pres- 
sure firora  the  ntolton  interior*  Water  sprays  surroimd  these  sizing 
rolls  and  the  ingot  is  quite  cool  irtien  it  gets  down  to  the  pinch  rolls. 
Williams  frequently  uses  mold  lubrication  to  increase  speed  of  cas- 
ting. The  lubricant  is  introduced  at  the  top  along  the  mold  wall, 
the  quantity  being  controlled  by  the  speed  of  the  pinch  rolls  drawing 
the  continuous  billet  out* 

The  Williams  machine  has  been  given  no  steady  trial  over  a  long 
run  on  steeli  however,  the  experimental  ingots  produced  have,  for  the 
most  part,  satisfactory  structure* 

A  number  of  makers  of  carbon  steel  ^o  habitually  handle  small 
heats  have  sliown  interest  in  the  process*  It  is  also  conceivable 
that  those  makers  of  high  alloy  steels  and  certain  types  of  stainless 
steels  may  see  in  this  process  a  possible  means  for  eliminating  present 
high  scrap  losses  and  a  more  flexible  method  of  handling  small  tonnage 
orders* 
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rETEEMdATIOH  OF  THE  TEMKHATUEE  AT  AHI  POIST  IM  THB  IHQOT 


In  the  following  investigation  certain  assumptions  are  made  to 
simplify  the  problem  and  bring  it  within  the  realm  of  possible  solution. 
These  assumptions  are  as  follows t 

(a)  Thenaal  conductivity  does  not  vary  with  tern-- 
perature» 

(b)  No  recognition  of  the  heat  of  fusion  of  steel 
is  oade  in  the  teinperature  distribution  in 
the  ingot* 

(c)  Conrection  in  the  molten  metal  is  neglected. 

It  is  realized  that  the  first  two  assunqptions  will  cause  appreciable 
error  in  the  final  results^  however,  unless  these  assa^tions  are  made 
the  problem  would  be  beyond  the  scope  of  this  paper.  The  assuaiption 
of  the  absence  of  convection  in  the  molten  steel  is  justified  since  it 
will  become  "raushy*  upon  entrance  into  the  mold  and  resist  any  ap- 
preciable turbulence. 

For  purposes  of  this  analysis  a  cylindrical  ingot  is  used  with 
all  results  in  cylindrical  coordinates.  Consider  first  an  elemental 
▼oluioe  as  follows  s 


The  partial  differential  eqoation  expressing  the  heat  oonduoti^ 
for  this  problem  is?  (refer  to  Appeadix  ^A*  for  derivation)  • 


uliero  o<   is  the  »tl»r«al  diffosivlty*'  of  U»  material* 

1!he  tl^naal  diffasd.'Tity  is  a  property  of  the  substance  in  iihidii 
heat  is  transferred  and  stored.  It  is  the  ratio  of  a  quantity  which 
is  proporti<Hial  to  the  heat  eondaotion  (numerator)  and  to  a  quantity 
wfaieh  is  proportional  to  the  he&t  storage  (denominator)*  lleat  unite 
eanoel  and  the  dimei^ons  become  ft^/hr* 

Aasoning  it  is  possible  to  express  the  soltstion  of  the  above 
equation  as  the  product  of  3  functions,  one  of  which  is  a  function  of 
•l*  alone,  one  a  function  of  ^   alone,  and  the  other  a  fonetion  of 
•a*  alone,  then 

The  9o3jatlon  (as  shown  Appendix  "B")  is 

latere  "a*  is  a  constant 

Ifdtei  Sisee  the  ten^ratiu^  distribution  Is  a  ftmoticm  of 
*»r"  and  *«*  only,  tenas  including  «d?*  drop  out* 
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fh»  Ingot  is  eonsidttred  to  be  a  sesrL-infinite  cylinder  id.th  the 
folloiring  boundary  conditions  s 

1.    ±-  iC<"^  ,^^^  .    o<r<R 

3.   t-^^^^Cz)  ,  o<«<-o,  r  =  ^ 

From  the  second  boimdary  condition  it  is  seen  that  for  t"?^  •** 
at7  =  *o  ,  C^  must  be  zero*  The  equation  then  reduces  to 

Solution  for  the  above  is  made  by  assuming  that  the  temperature  "t"  is 
the  sum  of  two  solutions,  that  is 

tj  «  t,  (z,r)         such  that  when  i-^^~  o  054<«»«,  rs  R 
t^=  t,  (2,r)         such  that    ft^'O  .  2=^  o  ,  os  r  <  K 

Solutions  for  t»  and  t^  are  taken  as  modified  forms  of  equations 

31  and  32 J  pp  192,  Carslaw  and  Jaeger,  Conduction  of  Heat  in  Solids  or 


where 


<,-  k-WM. 


and  the  a  ^  are  the  positive  roots  of  Jj!Ra^  s  0. 
The  total  solution  then  becomes  the  sura  of  equations  (1)  and  (2)  above. 
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Redacing  the  latter  part  of  equation  (2)  to  more  usable  form,  we  have. 


This  further  reciucea  to 


eHperature  distribution  in  the  ingot  thus  becomes 


For  actual  temperatures  the  functions  f(r)  and  l(/$   raust  be 
supplied  from  i*iysical  aspects  of  the  problem.  Consider,  for  steady 
flow,  the  foUoving  assumptions) 

f  (r)  Z    Tp   (pouring  teiaperature) 
F(  )  ■  1L(e>where  L  is  the  length  of  the  mold. 
Substituting  these  values  and  integrating  we  have 
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vlth  farther  slnpllflcatlon,  this  reduces  to 


4=^ 


Since  K^  is  negative  for  all  values  of  a^  this  vill  converge  very 
rapidly. 
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IMVE^lGkTKM  OF  MXIMUM  COOLIHG  Wtm  lUM  MSIES,  BEUTION  TO  HMT 
TH/IKSFSB  COEFFICIENT,  AND  EXPHESSIOH  FOR  MOLD  LENGTH 

For  this  phase  of  the  probl«a  consider  the  ma^itude  of  the 
heat  transferred  from  the  ingot  to  the  cooling  water  with  the  as- 
siu^tion  that  -Uie  velocity  of  the  ingot  is  such  as  to  cause  final 
freezing  at  the  instant  the  ingot  emerges  fron  the  boIgU     Since  the 
nass  rate  of  flow  of  the  netal  is  constant,  the  heat  transferred  in 
imit  time  for  a  mold  of  lengtii  L  is  equal  to 

Q  a  4(hin    -    ^omO     BWsec 

^ere      m  =•     —7—  V^  Y^  z    velocity,  ft/^r 

D    s     diameter  ingot 
^  Z    daisity  steel 


"^  hin  .    H  +  C^  Tp  Jp  «    P«»*^  ten^erature 

■^  f    ^  H     =    heat  of  fusion 

hout  «  C-,T  Cp  5    Specific  heat 

'  T     M    mean  temperature  of 

ingot  cross-section 
at  mold  exit 

Siibstituting  these  values  into  the  above  equation,  the  total  heat 
transferred  is 

Dividing  throu^i  by  tXDL,  we  have 
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Since  "q"  and  "V/  are  lndep«adent  tiiere  is  obtained  an  escpresslon 

for  the  length  of  the  mold  in  terms  of  the  heat  transferred. 


L-  Y  ^**^^^VT)1 


From  the  basic  law  of  heat  conduction  for  a  hollow  cylinder, 
the  heat  transferred  from  the  ingot  to  the  oooUng  water  is 

and  by  integration  within  the  lisdte  r^  and  T2 

where 

ri  s  radius  of  ingot 

r2  s  radius  of  mold,  water  side 

ti  M  arerage  surface  temperature  of  ingot  within 

mold 
t2  s  average  water  teaperature  in  the  mold 

The  tem  "U"  is  a  composite  or  "overall"  coefficient  of  heat  transfer 

and  is  defined  as  follows! 

'-    -L-  .  j5        4- 


U  '  -^w  ^  kw  ^  ^^ 


where 

h^  s  surface  coefficient  of  heat  transfer,  water 

side 
k^  2  -Uiermal  conductivity  of  mold 
h^^  s  surface  coefficient  of  heat  transfer,  ingot 

to  mold 
1  :  mold  wall  thickness 


from  this,  it  is  seen  that 
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8ab8tltuting  this  into  ttie  equation  of  heat  conduction  and  solving 
for  hj  ,  we  have 

For  liquids  with  aeynold»s  niaabere  exceeding  lO^CKX)  the  expressitm 
for  hy  i0  given  1:17  the  Colbum  equation 

(equation  Ud,  pp  I68,  McAdams,  Heat  Transmission) 


Applying  this  data  to  a  specific  problem,  consider  the  fol» 

lowing  datas 

Diameter  of  ingot  k  inches 

Pouring  tenqperature^  T  28OOOF 

Length  of  mold,  L  3  ft. 

Freezing  temperature  276OOF 

T  2560OF 

Heat  of  fusion  of  steel  6$  Btu/# 

Specific  heat  steel  0.1075  Btu/# 

Mold  thickness  0.5  inch 

kw/1  1650  B/hr.fb.F 
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frtm  this 

Q  s  11.1  X  10^  Btu/hr 

q  r  3.5U  X  ICP  BtuAr,ft2 

Using  this  value  of  "q"  in  the  equation  for  the  determination  of 
the  surface  coefficient  of  heat  transfer  (ingot  to  mold),  h^   can  be 
determined*  With  this  value  and  values  for  hy  and  kw/l,  an  "overall" 
coefficient  for  heat  transfer  for  varying  Reynold's  numbers  is 
calculated.  Curves  of  h^,  h^,  k^/l  and  U  versus  Reynolds  numbers 
are  shown  plate  1,   pp  iM. 

The  results  of  these  curves  snow  that  the  overall  coef- 
ficient of  heat  transfer  is  dependent  upon  the  surface  coefficisit 
between  the  ingot  and  the  mold,  and  for  this  particular  set  of  con- 
diticKis,  water  flow  rates  in  excess  of  Reynolds  number  $  x  10**  would 
be  iaq^ractical.  Such  determinations  can  be  made  for  varying  sets 
of  conditions  giving  a  family  of  similar  curves,  and  from  these  flow 
rates,  mold  length,  and  surface  coefficients  of  heat  transfer  can 
be  predicted*  These  results  are  not  conclusive  as  much  depends 
upon  reasonably  accurate  determination  of  temperatures  wil^in  the 
ingot*  It  does,  however,  show  that  maximum  cooling  water  flow  rates 
are  d«p«ident  upon  the  amount  of  heat  that  can  be  transferred  from 
the  ingot  to  the  mold,  i&ich  in  turn  is  dependent  upon  the  correspond- 
ing surface  coefficient* 
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COMCIJJSIONS 

!•  As  pointed  out  pareviously  the  solution  for  the  t^aperature 
distribution  is  onXj  approximate  because  the  heat  of  fusion  of  the 
metal  has  been  neglected  in  the  derivation.  The  solution,  however. 
Is  useful  in  heat  transfer  problems  for  detertainiBg  t«iperat(ires  in 
cylindrical  bodies  moving  with  constant  velocity,  and  also  in 
predicting  ten^rature  distributions  in  continuous  casting  problems. 
Solutions  for  different  velocities  can  be  modified  by  jK*oper  selection 
of  the  functions  f (r)  and  Fte)  which  can  be  fairly  accurately  pre- 
dicted from  experimental  data* 

2«  In  the  investigation  of  maximum  water  flow  rates,  and  their 
relation  to  surface  coefficients  of  heat  transfer  between  ingot  and 
mold,  it  is  shown  that  the  over-all  coefficient  of  heat  transfer  is, 
in  the  lirait,  a  function  of  the  surface  coefficient  between  the  billet 
and  the  mold.  From  the  basic  law  of  heat  conduction  and  for  a  given 
size  ingot,  it  is  concluded  that  to  increase  the  heat  transferred  some 
means  imst  be  provided  to  increase  this  surface  coefficient.  In  the 
Williams  process  mold  lubrication  is  used,  which,  in  addition  to  re- 
ducing friction  a  resistance,  increases  the  surface  coefficient  of 
heat  transfer  by  making  the  heat  transfer  across  the  space  more  nearly 
that  of  pure  conduction.  It  is  recommended  that  for  the  lubricant 
an  oil  of  high  thermal  conductivity  be  used  •  the  optimum  between 
the  lubricating  effect  and  the  increased  heat  transfer  coefficient  is 
left  for  future  investigation. 
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After  determination  of  the  maxiiaaia  heat  transfer  i^te,  the 
optisnam  lengUi  of  the  mold  can  be  determined  from  the  relation 


L-  '^'t>^-<.V^Vl 


It  is  to  be  noted  here,  hovever,  that  the  above  equation  is  based  on 
the  assumption  that  final  feezing  takes  place  just  as  the  ingot 
emerges  frota.  the  mold.  More  realistically^  this  final  freezing  point 
takes  place  some  distance  below  the  mold.  That  is,  the  laold  length 
need  be  only  sufficient  to  establish  a  "skin"  around  the  sBolten  metal 
which  is  strong  enough  to  prevent  rupture  either  from  excessive  load 
or  from  ferrostatic  pressure  from  the  molten  interior.  Investigation 
of  minimum  mold  length  is  recommended  for  further  work. 
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HECOMMENDATIOTS  fOU   JUHTHER  ¥(®K 

Since  the  techniqae  of  continuous  casting  is  complex,  and 
since  time  liadtations  precluded  aimljrsis  of  all  pdiasef ,  the  following 
recommendations  are  inade  for  considerations 

!•  Frictional  cJrag  exerted  on  the  mold  hj  the  ingot 

Although  it  is  claiined  that  the  molten  metal  instantaneously 
freezes  upon  entry  into  the  mold  and  shrinks  away  from  frictional 
resistance  (The  Iron  Age,  April  h  and  April  11,  191*0),  it  seems  highly 
iiTiprobable*  One  of  the  lavs  of  sliding  friction  is  that  the  friction 
is  proportional  to  the  force  pressing  the  surfaces  together.  This 
force  is  increased  by  the  fcrrostatic  pressure  ftrom  the  molten 
interior,  "Hhich,  it  is  believed,  vill  overcome  the  shrinking  due 
to  freezing  of  the  nratal.  Adding  a  lubricant  to  decx^ease  the  coef- 
ficient of  friction  will  increase  the  maxiraim  load  capacity  and  hence 
the  maxiiBwa  allovable  velocity* 

Investigatiwi  of  the  aeehanlcal  friction  is  sin^lified  be- 
cause, within  limits,  sliding  friction  is  independent  of  the  area  of 
contact,  and  (except  at  start)  of  the  speed  of  relative  motion. 
The  problem,  therefor*,  resolves  itself  to  a  balance  of  forces  be- 
tween the  ferrostatie  pressure  and  the  ahrinklng  due  to  cooling. 
2.  MaxLiBom  Velocity  of  ingot  travel. 

The  velocity  of  the  ingot  vmries  directly  as  the  amount  of 
heat  which  can  be  satisfactorily  removed  from  the  steel.  With  suf- 
ficient heat  removal  flpom  the  mold  a  skin  will  be  developed  such  that 
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tlie  ingot  can  withstand  the  load  inposed  by  the  pinch  roUjj.  The 
maxiauia  velocity  will  be  determined  by  the  maximiiBi  heat  reaaoval 
consistent  with  strength  characteristics  of  the  skin  of  the  billet» 

3*   Load  capacity  of  ingot  and  determination  of 

ingot  breaking  point 

The  load  capacity  of  the  ingot  is  a  function  of  the  load 

capacity  of  the  skin  formed  in  the  nold.  The  load  will  be  caused 

by  surface  tension  and  surface  film  effects  of  the  steel  and  friction- 

al  resistance  to  the  passage  of  the  ingot.  This  load  c«gpacity  will 

reach  its  maximum  utien  such  velocity  is  reached  that  the  maxiiaum 

heat  removal  will  be  below  that  necessary  to  develop  sufficient  skin 

thickness  and  strength.  Since  the  ferrostatic  pre8e\ire  increases 

down  throu^  the  mold,  the  frictiooal  resistance  will  also  increase, 

and  hence  some  point  will  be  reached  where  the  load  exceeds  the 

ultimate  strength  of  the  billet.  Experimental  evidence  has  shown 

that  the  breaking  point  will  occur  some  distance  down  the  mold,  but 

the   exact  location  has  not  been  determined. 

U»   Investigation  of  possible  surface  boiling  on  the  water 
side  of  the  mold 

With  very  hi^  rates  of  heat  transfer  and  hi^  tenqjerature 

gradients  produced  in  the  mold  wall  there  are  excellent  possibilities 

of  surface  boiling  in  the  cooling  water.  With  hi^  feed  rates,  a 

small  fraction  is  vaporized,  and  the  boiling-side  coefficients  may 

be  roughly  twice  those  predicted  from  straight  water  cooling.  These 

increased  coefficients  may  be  due  to  increased  turbulence  in  the  film 

near  the  wall.  Investigation  along  ti:ieee  lines  can  be  very  instru- 
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mental  in  Increasing  heat  transfer  rates. 

5*   Surface  Tension  and  Surface  Flla  Sffects 
k  higb  value  of  surface  tension  has  the  effect  of  increasing 
the  pressure  required  to  cause  tae  metal  to  enter  and  flow  almig  a 
narrow  channel;  but  the  direct  effect  of  the  true  surface  tension 
is  alaaost  overshadowed  ly  the  effect  of  surface  films,  the  influetice 
of  which  is  equivaleeut  to  Va«t  of  a  greatly  increased  surface  tension* 
It  has  bean  found  ^lat  surface  OKlda  filas  af feet  the  flotfing  power 
of  metals  with  the  aj^^rent  surface  tension  of  some  metals  increased 
three  times*  Steel  has  a  very  high  value  of  surface  tension. 
(Brlgg8«  Metallurgjr  of  Ste«l  Castixigs) 

Boto  effects  will  add  to  the  resistance  to  flow  and  decrease 
the  max! mam  load  capacity  of  the  ingot.  Althou^  it  is  thou^t  to 
be  considerably  small  coo^>ared  to  frictional  resistance,  further 
investigation  ^euld  be  initiated  to  determine  its  effects* 
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APPENDIX    A 
gqufcticgi  of  Condactlon  for  a  MorLng  Cylinder 
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Beat  condttctcd  throiigh  "r"  face  is 


-d€?.. 


Heat  conducted  through  the  "r  -♦•  dr"  face  is 
Heat  conducted  through  "(J?"  face 


dQ,^ 


Note J  Temperature  gradient  in cpdirection  is  the 
diange  of  temperature  with  respect  to  ihe 
distance  represented  hy  the  angle^^  whidi 
distance  also  depends  on  the  radius  r. 


Heat  COTiducted  through  "   cP  +  q/  cP    "  face 
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^^'^^h)^' 
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Heat  conducted  througa  "a"  face 

^^., 

r  v?^ 


-^f-s.^ 


Heat  cwiducted  tiirou^  "z  +dz"  face 

or      k^  '-c/^^<^'C+  kd<Cdrdir^,dr-  •   --   d(5^^ 

Heat  gained  liy  ihe  elemental  "roliaae  due  to  caiduction  is  equal  to 

ZdQ^~    ^dc3, 

iirtiich  equals 

Heat  stored  in  tne  dif feraitial  body 

However  the  temperature  incremaat  of  a  particle  of  trie  fluid  in 
time  interval  dT  becomes  .  i 
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&iiiQ«  «•  are  considering  veleolty  in  the  8*»dir6etlaci  and  a  ecNadLtl<»i 
of  0t«ad^  state*  this  reduces  to 

Ftm  thi8»  tlieUf  the  Heat  stored  la  the  bo<!^  becomes 

Since  fio  heat  snergsr  is  being  developed  in  the  boc^*  ttM  lAtf.of  the 
eonserfation  of  energy  reduces  to 

lieat  stored  s  heat  in  •>  heat  out 

Substituting  and  dividing  throu|^  by  rcl<0^  <^^2^aT' 

iAtrodofiiiM^  the  t«m  "theraal  dif fosivitT* 

(fhe  thennaX  dif ftisivit/  is  a  property  of  the  eubstanoe  in  which 
heat  is  iranitf erred  and  stored)* 

Sxibetituting  this  value >  the  equation  1 
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APPQJDIX  B 
SOLUTION  FOR  THE  SQUATIOM  OF  CQKDUCTIQH  »  TKMP^iiiATlIRE  DISTRIBUTION  IN 
INGOT 

AsBume  it  is  possible  to  express  the  solution  of  the  dif* 
ferential  equation  (Appendix  A)  as  the  product  of  three  functions « 
one  of  which  is  a  functicai  of  "r"  alone,  one  a  function  of  "<^« 
alone,  and  the  other  a  function  of  "z"  alone,  then 

The   equation  takes  the  fom, 

dividing  throu^i  bj  R  J)-?*  *  transposing  and  equating  to  a  constant, 
«•  hav« 

FlroB  the  left  side  of  the  equati<m,  vrlte 

Which  gives 

the  solution  of  which  ±9 
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Bquation  (2)  also  gives 

lot  X  :  ar 

tiieaa  since  , 

eqiiation  (^)  becomes 

This  equaticxi  is  Be8sel*e  diffsrential  eqiiati<Hi  (Hildebrands  Ad-vanced 
Calculus  of  Engineers,  pp  l65»l67)>  and  the  genei^l  soluticsn  for  "p" 
an  integer  can  be  written 

tAiere  Ep  (x)  is  giv^i  as  follovs  (Woods,  Advanced  Calculus,  pp  288) 

n-o  •  h-o  ^    ' 

Since  equation  (?)  is  applied  to  a  solid  cylinder,  Bessel^s  function 
Kp(x)  becomes  infinite  at  x  r  0  (in  0  s~*>o)^  and  it  follows  that 
C2  2  0.     Thus  for  the  solid  cylinder,  the  general  solution  is  taken 
asi 
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Taking  the  remaining  equation  in  otir  solution 

Using  the  operator  method  of  solution,  the  equation  becomes 

[t5"-l-D-a]z.-  o  do) 

The  solution,  therefore,  of  equation  (9)  becomes 
uhere 


-=  h  ^ 

\tH 

2.        o 

Uhi/  +  0^ 

0=^ 

Haace  for  the  cylinder  the  solution  of  the  differential  equation  is 
given  by  substituting  equations  (ii.),  (8),  and  (ll)  into  equation 

Since  the  temperature  distribution  is  a  function  of  "r"  and  "z"  only, 
thffli,  p  «  0,  and  the  solution  of  the  equation  becomes 
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